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We have performed for the first time a comprehensive study of the sensitivity of r-process nu¬ 
cleosynthesis to individual nuclear masses across the chart of nuclides. Using the latest version 
(2012) of the Finite-Range Droplet Model, we consider mass variations of ±0.5 MeV and propagate 
each mass change to all affected quantities, including Q-values, reaction rates, and branching ratios. 

We find such mass variations can result in up to an order of magnitude local change in the final 
abundance pattern produced in an r-process simulation. We identify key nuclei whose masses have 
a substantial impact on abundance predictions for hot, cold, and neutron star merger r-process 
scenarios and could be measured at future radioactive beam facilities. 

PACS numbers: 21.10.Dr,26.30.Hj,25.20.-x,26.30.-k 


I. INTRODUCTION 

One of the most challenging open questions in all of 
physics is the identification of the site or sites of rapid 
neutron capture, or r-process, nucleosynthesis mm- Pro¬ 
duction of the heaviest r-process elements requires on the 
order of 100 neutron captures per seed nucleus; exactly 
where and how such rapid neutron captures occur has 
yet to be definitively determined [3]. 

One attractive potential site is within the cold or 
mildly heated tidal ejecta from neutron star or neutron 
star-black hole mergers [1] . Current state-of-the-art sim¬ 
ulations show a vigorous r process with fission recycling 
in the merger ejecta [SHE]. The resulting abundance 
pattern is relatively insensitive to variations in the ini¬ 
tial conditions, which naturally explains the consistent 
56 < ^ < 82 pattern observed in the solar system and 
r-process-enhanced halo stars liiin]. Vigorous produc¬ 
tion of radioactive r-process nuclei can also lead to an 
observable electromagnetic transient accompanying the 
merger event [IIHI3], an example of which may have al¬ 
ready been detected m [n]. It is less clear whether 
mergers happen often enough or early enough in galactic 
history to fit all of the observational data [TH1 - IT3] . The 
neutrino-driven wind within a core-collapse supernova is 
perhaps the best-studied alternative HHHD], though the 
combination of moderate neutron-richness, high entropy, 
and fast outflow timescale required to make the heavi¬ 
est r-process elements does not appear to be achieved 
in modern simulations Instead attention has 

shifted to more exotic sites connected to the deaths of 
massive stars, including neutron-rich jets [25], supernova 
neutrino-induced nucleosynthesis in the helium shell |26j , 
and collapsar outflows EZlllH]- 

In principle, the proposed environments have such dis¬ 
tinct astrophysical conditions that their abundance pat¬ 
tern predictions should be clearly distinguishable. Cur¬ 
rently simulations lack this precision |29l[30] , in large part 
due to uncertainties in the required nuclear data. Proper¬ 


ties such as masses, neutron capture rates, /3-decay rates, 
and /3-delayed neutron emission probabilities are needed 
for thousands of neutron-rich nuclear species from the 
valley of stability to the neutron drip line. Presently 
there is little experimental information available for the 
vast majority of these quantities. Simulations must in¬ 
stead rely on extrapolated or theoretical values, where 
different approaches can produce markedly different (and 
often divergent) predictions. 

Nuclear masses are particularly important for the r 
process as they enter into the calculations of all of the 
aforementioned nuclear properties, which shape how each 
phase of the r process proceeds. In a classic r process, 
an equilibrium is established between neutron captures 
and photodissociations, and nuclear masses directly de¬ 
termine the r-process path through a Saha equation: 

Y{Z,N + l) ^_G{Z,N + l) Nr, __ lSr,{Z,N+iy 

Y{Z,N) 2G{Z,N) (fcT)3/2®^P[ fcT 

( 1 ) 

where G{Z, N) are the partition functions, Nn is the neu¬ 
tron number density, kT is the temperature in MeV, and 
Sn{Z,N + 1) is the neutron separation energy, the differ¬ 
ence in binding energy between the nuclei {Z, A^± 1) and 
{Z, N). Each isotopic chain is connected to its neighbors 
by /3-decay, and thus the /3-decay lifetimes of nuclei along 
the r-process path set their relative abundances. Modern 
nuclear network calculations show that this equilibrium 
picture is an excellent approximation for early-time r- 
process evolution in many astrophysical scenarios. Even¬ 
tually (n, 7 ) ^ ( 7 , n) equilibrium fails, or in some sce¬ 
narios is not established at all, and then neutron capture, 
photodissociation, and /3-decay all compete to shape the 
final abundance pattern. 

The roles of these individual pieces of data in r-process 
nuclear network simulations have been examined via sen¬ 
sitivity studies. In these studies, baseline astrophysical 
conditions are chosen, a single nuclear property is var¬ 
ied, and the simulation is rerun with the nuclear data 
change and compared to the baseline. Sensitivity stud- 
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ies highlight the pieces of data with the most leverage 
on the final abundance pattern and elucidate the mecha¬ 
nisms of influence. They have so far been performed for 
neutron capture rates , /3-decay rates [Ml |M] , and 

photodissociation rates via their dependence on nuclear 
masses [3M38] . and are in progress for /3-delayed neu¬ 
tron emission probabilities |39j . These studies are dis¬ 
tinct from and complementary to studies of the influence 
of groups of nuclear properties, e.g. gniiii], global the¬ 
oretical model predictions, e.g. [42], or global variations 
in a Monte Carlo approach [521 [30] that can quantify 
correlations between nuclear physics inputs |43| . 

The sensitivity studies described above looked for the 
individual nuclear properties with the greatest impact 
on the r process by considering variations of one piece 
of data at a time. However, we know that modifying 
a single nuclear mass alters all of the nuclear proper¬ 
ties that depend on that mass. Here we perform nuclear 
mass r-process sensitivity studies in which variations in 
individual nuclear masses are consistently propagated to 
all affected nuclear properties, including neutron capture 
rates, photodissociation rates, and /3-decay properties. 
This approach was developed in Ref. [M] for spherical nu¬ 
clei and is extended for the first time in this work to the 
entire chart of the nuclides between Z = 30 and Z = 80. 
These new studies capture the full impact of the uncer¬ 
tainties in individual masses on the r process and point 
to the most important masses to measure in present and 
future experimental campaigns. Additionally they place 
direct constraints on the precision needed for measure¬ 
ments of nuclear properties and theoretical models, in 
order to improve r-process predictions and, eventually, 
distinguish between possible astrophysical sites. 


II. BASELINE r-PROCESS SIMULATIONS 

We begin our sensitivity studies with a choice of base¬ 
line astrophysical trajectories. We investigate a range 
of astrophysical conditions including a low entropy hot 
wind, high entropy hot wind, a cold wind and a neu¬ 
tron star merger. In a hot wind an equilibrium is es¬ 
tablished between neutron captures and their inverse re¬ 
action, photodissociation, and the /3-decays that move 
the nuclear flow to higher atomic number, Z, control the 
timescale for heavy-element production. Nuclear masses 
are influential during equilibrium as they directly set the 
r-process path for a given temperature and density, as 
shown in Eqn. To explore nuclear mass uncertainties 
in a low entropy hot wind we use a parameterized wind 
model from |45| which has an entropy of 30 fcs, an elec¬ 
tron fraction of Ye = 0.20 and a timescale of 70 ms. This 
trajectory yields the production of heavy elements out to 
the third peak but is not neutron rich enough to lead to 
fission recycling. The high entropy hot trajectory uses 
the same parameterization with entropy of 100 ks, an 
electron fraction of Yg = 0.25 and a timescale of 80 ms. 
In a cold wind the (n, 7 ) ^ ( 7 , n) equilibrium is short¬ 


lived as the temperature drops quickly and photodisso¬ 
ciation becomes negligible. The r-process path moves 
far from stability where a new quasi-equilibrium can be 
established between neutron captures and /3-decays. To 
study nuclear masses under these conditions we employ 
a neutrino-driven wind simulation of Ref. [46] with ar¬ 
tificially reduced electron fraction of Ye = 0.31 to pro¬ 
duce a main r process. We also consider a neutron star 
merger environment using a trajectory from Bauswain 
and Janka, for which the mildly-heated ejecta is suffi¬ 
ciently neutron-rich to undergo fission recycling [^ . 

In all cases, once the supply of free neutrons is con¬ 
sumed the r-process path begins to move back to sta¬ 
bility. The criterion of neutron exhaustion signals the 
start of the freeze-out phase of the r process in which 
key abundance features are formed, such as the rare earth 
peak [47l [48]. Additional neutrons during this time come 
from photodissociation, neutrons emitted promptly after 
/3-decay, or fission [5^. During this freeze-out phase the 
reaction flows are sensitive to individual masses of nuclei 
between the r-process path and stability. 

Our studies employ a dedicated r-process reaction net¬ 
work code gziisni which has been updated in recent stud¬ 
ies [laiM] and includes a schematic treatment of fission 

[52] . For the baseline nuclear masses we use the 2012 
version of the Finite-Range Droplet Model (FRDM2012) 

[53] . This model has an root-mean-square (rms) error of 
0.57 MeV compared to known masses from the Atomic 
Mass Evaluation (AME2012) [54]. When available, we 
use measured masses from the AME and keep these 
masses unchanged during our calculations. The use of 
experimental (AME2012) and theoretical (FRDM2012) 
mass values requires some care to avoid large disconti¬ 
nuities. Therefore when calculating quantities that de¬ 
pend on two masses of different origin (theory and ex¬ 
periment), we make sure to consistently calculate that 
quantity within a given dataset. For example, if a mea¬ 
sured mass is available for {Z + 1,N — 1) but not {Z, N), 
the /3-decay Q-value for {Z, N) is calculated from the 
FRDM2012 masses only. 

The baseline neutron capture rates are calculated with 
the publicly available statistical model code TALYS [55] 
using FRDM2012 and the latest compilation of mea¬ 
sured masses from the AME2012 as described above. 
We use the default settings for level density, 7 -strength 
function and particle optical model. The default level 
density model used in TALYS combines a constant nu¬ 
clear temperature at low energies and matches it with 
a back-shifted Fermi gas model using systematics [56] . 
Masses can enter into this level density model via the 
definition of the constant temperature, which, as used 
in the TALYS code, is proportional to one over the 
square root of the shell correction term: dW{Z, N) = 
M{Z,N) — Mlom{Z, N) where M{Z,N) is the mass of 
the nucleus and Mldm{Z, N) is the predicted mass to 
a spherical liquid-drop. Additionally, the level density 
parameter, a, is proportional to dW and hence the nu¬ 
clear masses. The default 7 -strength function used in 
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TALYS is the formulation set out in Kopecky-Uhl (KU) 
m- Nuclear masses enter into this parameterization of 
the giant dipole resonance via nuclear temperature term 
which prevents the 7 strength from going to zero as 7 
energy decreases. The nuclear temperature itself is pro¬ 
portional to the square root of the neutron separation 
energy, S'„, as well as a, evaluated at S'„. 

Photodissociation rates are calculated from neutron 
capture rates by detailed balance: 


X-y{Z, N) oc exp 


Sn{Z,N)l 

kT 


(^z,N- 1 ) 


( 2 ) 


where Sn{Z, N) = M{Z, N—1)—M{Z, N)+Mn is the one 
neutron separation energy, M{Z, N) and M{Z, N—1) are 
masses of the nuclides, M„ is the mass of the neutron, T is 
the temperature, {crv) {z,N-i) is the neutron capture rate 
of the neighboring nucleus and k is Boltzmann’s constant. 
The baseline photodissociation rates use FRDM2012 & 
AME2012 masses and the neutron capture rates as cal¬ 
culated above. 

For the baseline ,5-decay rates we use experimental val¬ 
ues [SS] where available and theoretical estimates every¬ 
where else. Theoretical rates are calculated as in [59]: 




t 


1/2 


(3) 


where Xjs is the /3-decay rate, ti /2 is the half-life, i de¬ 
notes the ith excited state of the daughter nucleus with 
energy Ei, uji = (Qp — Ei)/'me is the /3-decay energy 
to this state in units of electron mass, I is the type of 
the decay, either Gamow-Teller (GT) or First-Forbidden 
(FF), / is the phase space factor and C is the /3-strength 
function. We use the /3-decay strength data from |5U] 
and calculate the phase space piece using FRDM2012 & 
AME2012 masses as described above. Neutron emission 
probabilities are calculated by combining the Quasipar¬ 
ticle Random-Phase Approximation (QRPA) model of 
Ref. |59| with a Hauser-Feshbach (HF) model from Ref. 
[ 60 j which has been extended to the neutron dripline in 
this work. This combined QRPA-HF approach produces 
larger predictions of average neutron emission towards 
the dripline compared to older (QRPA) methods, as it 
follows the statistical decay until the initial available ex¬ 
citation energy is exhausted. 


III. MASS VARIATION PROPAGATION 

For the sensitivity studies of this work we consider 
mass variations of ±0.5 MeV, approximately equal to the 
rms error of FRDM2012. Since measured masses tend to 
have experimental uncertainties much smaller than this, 
we restrict our individual mass variations to nuclei with 
extrapolated or unknown masses. 

When an uncertain nuclear mass is varied we recal¬ 
culate all the relevant nuclear properties of neighboring 
nuclei that depend on the changed mass. Specifically, if 
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FIG. 1: Shows the quantities of neighboring nuclei of im¬ 
portance to the r-process that may be altered by a change in 
mass of nucleus with Z protons and N neutrons. 


the mass of a nucleus {Z,N) with Z protons and N neu¬ 
trons is varied then it can lead to changes in the neutron 
capture rates of {Z,N) and {Z,N — 1), the separation 
energies of {Z,N) and {Z,N ± 1), the /3-decay rates of 
{Z,N) and {Z — 1,N + 1), and /3-delayed neutron emis¬ 
sion probabilities oi{Z,N), {Z —1,N + 1), {Z —1,N + 2), 
up to {Z — 1,N + 12) as shown in Fig. 

To propagate the changes to the neighboring neutron 
capture rates we invoke the ‘massnucleus’ command in 
TALYS and continue to use the default settings for level 
density, 7 -strength function and particle optical model 
with the varied mass. A ±0.5 MeV mass variation re¬ 
sults in a change in neutron capture rates of a factor of 
approximately two to five. Note that this corresponds to 
the uncertainty in the neutron capture rates due only to 
uncertain nuclear masses; larger variations in the rates 
come from choice of statistical model components and 
the treatment of direct capture m- The impact on the 
r process of these larger neutron capture rate uncertain¬ 
ties are addressed in earlier studies |32l |33l |50| . 

For the /3-decay rates, all measured half-lives remain 
unchanged during our mass variations and we adjust only 
the affected theoretical rates. For the theoretical values, 
the majority of the dependence on nuclear mass is con¬ 
tained in the phase space piece, which goes as /(w) ~ w® 
for allowed decays. This is a significantly stronger de¬ 
pendence on masses than appears in the nuclear matrix 
elements calculated in the QRPA [33]. Therefore when 
updating the theoretical rates with the mass variations 
we recalculate only the phase space factors and leave the 
/3-strength functions unchanged. This greatly reduces 
the numerical cost of the calculation while still capturing 
most of the dependence of the half-lives on the masses. 
We find mass variations of ±0.5 MeV lead to changes in 
the half-lives by a factor of roughly two to four. 

We recalculate the affected /3-delayed neutron branch¬ 
ing ratios in our QRPA-HF approach, assuming the 
structure of the ground state does not change with the 
mass variation. Modifications come from a change in the 
/3-decay Q-value of the parent nucleus or from a change 
to the neutron separation energies of the daughter nuclei. 
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A general dependence of neutron emission probabilities 
on changes in mass using QRPA-HF is thus entangled 
with the inclusion of 7 -ray competition at each neutron 
emission stage. 

IV. MASS SENSITIVITY STUDY 

With the required baseline nuclear data generated as 
described in Sec. [IT] and the computational tools in place 
to propagate the mass variations to the affected nuclear 
properties as described in Sec. |III[ we begin our sensi¬ 
tivity studies with the baseline astrophysical trajectories 
identified in Sec. [IT] For each nuclear species {Z,N) we re¬ 
peat the baseline r-process simulation twice: once where 
the mass M{Z,N) is increased by 0.5 MeV and once 
where the mass is decreased by 0.5 MeV. We then com¬ 
pute the metric 

F = 1Q0J2\X{A)-X,{A)\ (4) 

/I 

where Ab(A) is the final isobaric mass fraction in the 
baseline simulation, X (A) is the final isobaric mass frac¬ 
tion of the simulation when all nuclear inputs have been 
modified based off the change in a single mass, and the 
summation runs over the entire baseline pattern |44j . As 
defined here, the F metric quantifies the impact of the 
mass on the global r-process pattern. This procedure is 
then repeated for every nuclear species from the limits 
of the AME2012 experimental values to the FRDM2012 
neutron drip line. 

Fig. [2] shows the results of the full-chart mass sensi¬ 
tivity studies run for our four choices of astrophysical 
conditions. The shading of each rectangle represents the 
largest of the two F values that we obtained when we cal¬ 
culated the r-process abundances for the reference mass 
(FRDM2012) (I) plus 0.5 MeV and (2) minus 0.5 MeV. 
Roughly 1100 nuclei are included in each study with 
30 < .Z < 80 and 60 < A < 130. As a general result, we 
find influential nuclear masses lie along the equilibrium 
r-process path as well as along the decay pathways back 
to stability, most of which are within the predicted exper¬ 
imental reach of the upcoming Facility for Rare Isotope 
Beams (FRIB) [62]. This is particularly evident around 
the closed neutron shells where we find the largest global 
impact of nuclear masses. 

Qualitatively, the overall distribution of influential nu¬ 
clear masses in the four cases reinforces the conclusions 
from our previous work [MlEHlIllI- In particular the pat¬ 
tern of most impactful nuclei for the hot wind r process 
shown in the top panel of Fig. is directly compara¬ 
ble to that from [36l |38] . Our new studies show several 
key improvements, which come from addressing the lim¬ 
itations of these early studies. Refs. jSS] I3E| considered 
the propagation of mass variations to photodissociation 
rates only, which underestimates the resulting sensitiv¬ 
ity measures, especially in cold wind or merger r-process 
scenarios where photodissociation is suppressed. In these 


cases, the influence of masses on the r process occurs via 
the decay properties and neutron capture rates—an ef¬ 
fect which was captured for the first time in Ref. jUj but 
only for closed shell nuclei. Thus our new studies are 
the first reliable full-chart estimates of the impact of nu¬ 
clear masses in cold and merger r-process scenarios. In 
addition, the hot wind studies show higher mass sensi¬ 
tivities in the rare earth region compared to the earlier 
studies, due to the interplay of /3-decay, neutron cap¬ 
ture, and photodissociation that forms the rare earth 
peak sziiini. The exact mechanisms by which each piece 
of nuclear data influences the final abundance pattern in 
different types of r processes are explored in Refs. ini- 
ESIElllMlEHlISlllil]. 

Any nucleus with impact parameter above A ~ 10 
(purple shading) implies a significant global change to the 
abundances. Only Palladium ( 46 Pd), Cadmium ( 48 Cd), 
Indium ( 49 ln), and Tin (soSn) isotopes show sensitivities 
above A = 20 in our studies. Table 1 includes the F 
measures of these most impactful nuclei. 

TABLE I: Important nuclei from Fig. with Fmax > 20. 


z 

N 

A 

P 

■*- max 

Trajectory 

48 

84 

132 

101.39 

high entropy hot 

50 

86 

136 

83.68 

high entropy hot 

49 

84 

133 

74.59 

high entropy hot 

49 

86 

135 

73.82 

high entropy hot 

49 

85 

134 

72.84 

high entropy hot 

49 

87 

136 

71.04 

high entropy hot 

49 

88 

137 

68.08 

high entropy hot 

49 

89 

138 

66.43 

high entropy hot 

48 

88 

136 

58.42 

low entropy hot 

46 

86 

132 

55.56 

nsm 

48 

86 

134 

52.38 

low entropy hot 

46 

86 

132 

48.67 

cold 

46 

84 

130 

43.73 

cold 

48 

90 

138 

37.86 

low entropy hot 

46 

85 

131 

34.81 

cold 

50 

88 

138 

29.64 

high entropy hot 

46 

84 

130 

24.19 

nsm 

48 

88 

136 

24.10 

cold 

48 

85 

133 

23.38 

low entropy hot 

48 

87 

135 

21.12 

low entropy hot 

48 

89 

137 

20.63 

cold 


V. ABUNDANCE PATTERN PREDICTIONS 

We now use the above sensitivity study results to ex¬ 
amine the variances that can arise in the final abun¬ 
dance pattern due to uncertainties in nuclear masses. 
In each sensitivity study we first select the nuclei for 
which F > 0.1 as shown in Fig. We create an en¬ 
semble of abundance patterns from the simulations that 
include the individual mass changes ±0.5 MeV of these 
nuclei. We then compute the variance of this reduced en¬ 
semble of abundances for each value of A. The variance 
in the final abundance patterns for the low entropy hot 
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FIG. 2: Nuclei that significantly impact final r-process abundances for four astrophysical conditions. The intensity of the color 
denotes the maximum F value, Eqn. resulting from individual ±0.5 MeV mass variations. Light gray denotes the extent 
of measured masses from the 2012 AME and stable nuclei are colored black. For reference, estimated accessibility limits are 
shown for the upcoming Facility for Rare Isotope Beams (FRIB) (black line - intensity of 10“^ particles per second [63]). 
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FIG. 3: Variances of the ensembles of final abundance pat¬ 
terns (shaded bands) for the four sensitivity studies described 
in this work, compared to scaled solar r-process residuals 
(black circles) from |64| . 


wind, high entropy hot wind, cold wind and neutron star 
merger studies are shown in Fig. [^compared to the solar 
isotopic r-process residuals. 

The dependence of r-process predictions on the uncer¬ 
tainties in nuclear masses is shown by variance bands for 
four different astrophysical conditions in Fig. These 
four trajectories have distinct r-process paths and dy¬ 
namics during freeze-out which means the variance in 
abundances in each environment comes from a different 
aspect of the dependence on nuclear masses as mentioned 
above. This result is only obtainable by using our ap¬ 
proach of consistently propagating uncertainties from nu¬ 
clear masses to all of the relevant nuclear quantities for 
the r process. 

One may be tempted to rule out the hot conditions 
as the variance bands are clearly offset for the third 
{A = 195) peak. However, we note that our approach 
here actually underestimates these bands in all scenar¬ 
ios — larger uncertainties from reaction rate calculations 
are not included, which are particularly important for 
the hot scenarios, and mass changes in these studies are 
done on an individual basis. Methods that rely on a 
global Monte Carlo approach have the ability to resolve 
these drawbacks. Preliminary work |29j in this direction 
suggests details of the abundance pattern can be clearly 
resolved if mass uncertainties are reduced to less than 
0.1 MeV. The next generation of nuclear mass measure¬ 
ment campaigns will be crucial in the progress toward 
this ambitious goal. 

With this caution in mind, Fig. [^confirms our previous 
results that mass model uncertainties are currently too 
large for precision abundance pattern predictions capa¬ 


ble of differentiating between r-process conditions |44) . 
For example, the formation of the A ~ 160 rare earth 
peak can in principle be used to constrain the r-process 
site [3H]j however the variance bands in Fig. [^are larger 
than the peak itself. This indicates that the features of 
the mass surface in this region responsible for rare earth 
peak formation are likely on the order of the rms value 
of FRDM2012 or smaller. 

We also note that the new FRDM2012 masses show 
marked improvement over the FRDM1995 masses in 
matching features from the solar pattern, as is clear from 
a comparison of the final abundances of Fig. to those 
in Ref. [U]. The most notable improvement of the mass 
model comes from an enhanced description of nuclei in 
the transition region between the N ^ 82 region and the 
rare earth region. 


VI. CONCLUSION 

In summary, we have shown for the first time how un¬ 
certainties in individual nuclear masses propagate to in¬ 
fluence and shape the r-process abundance distribution 
across the chart of nuclides. We consider variations of 
individual nuclear masses and recalculate consistently all 
relevant Q-values, neutron capture rates, photodissoci¬ 
ation rates, /?-decay rates and /3-delayed neutron emis¬ 
sion probabilities, as shown in Fig. We find mass 
uncertainties of ±0.5 MeV have a significant impact on 
r-process abundance predictions as summarized in Fig.[^ 
In terms of our metric, a value of F ~ 20 represents a 
large local change or a global shift in r-process abun¬ 
dances. 

We explore changes to masses of ±0.5 MeV from 
FRDM2012 in four astrophysical trajectories: a low en¬ 
tropy hot wind, a high entropy hot wind, a cold wind, 
and a neutron star merger. Shifts in the equilibrium 
path play the dominant role in a hot r process, and our 
results here mirror earlier studies [36l [38] where mass 
variations were propagated only to the photodissocia¬ 
tion rates. Changes to weak decay properties and neu¬ 
tron capture rates are essential to include particularly in 
the cold wind and merger cases where photodissociation 
channel is suppressed. We find a similar dependence on 
nuclear masses for all astrophysical conditions studied, 
as shown in Fig. due to the propagation of mass un¬ 
certainties to all relevant quantities. The nuclei with the 
most impactful masses lie along the equilibrium r-process 
path, as expected, and also along the decay paths to 
stability. This result strongly reinforces our conclusions 
from previous studies that understanding freeze-out is 
critical for predicting r-process abundances under any 
astrophysical conditions. 

Many of the influential nuclear masses identified in our 
studies will be accessible to future radioactive beam facil¬ 
ities. These measurements have the potential to dramati¬ 
cally improve the precision of r-process simulations. Still, 
some masses will remain beyond experimental reach. It is 
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therefore also of key importance to improve global mass 
models. This work provides additional motivation for 
strengthening efforts in improving the global description 
of nuclear masses as it clearly shows that current mass 
model uncertainties are still too large for the finer de¬ 
tails of the isotopic r-process abundance pattern to be 
resolved. The question then arises: how well do we need 
to know nuclear masses in order to predict a final com¬ 
position of the r process? Studies which address this 
question strongly point to the resolution of abundance 
features if global model uncertainties are reduced to less 
than 0.1 MeV [5S1 Thus, a concerted effort is re¬ 
quired by the community to analyze model inputs and 


physical assumptions in order to achieve reliable calcula¬ 
tions as close to this threshold as possible [65H67] . 
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